Magnetic order in ferromagnetically coupled spin ladders 
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A model of coupled antiferromagnetic spin- 1/2 Heisenberg ladders is studied with numerical tech- 
niques. In the case of ferromagnetic interladder coupling we find that the dynamic and static 
structure factor has a peak at {n,n/2) where the first (second) direction is along (transversal) to 
the ladders. Besides, we suggest that the intensity of this peak and the spin-spin correlation at 
the maximum distance along the ladder direction remain finite in the bulk limit for strong enough 
interladder coupling. We discuss the relevance of these results for magnetic compounds containing 
ladders coupled in a trellis lattice and for the stripe scenario in high-Tc superconducting cuprates. 

PACS: 75.10.Jin, 75.40.Mg, 74.72.-h, 75.50.Ee 



I. INTRODUCTION 

One of the main topics in condensed matter physics 
in recent years has been the study of low-dimensional 
antiferromagnetic spin systems. The strong interest in 
this field has been sparkled by the realization that Cu02 
planes play an essential role in high-Tc superconductors, 
which was followed by the appearance of several com- 
pounds characterized by the presence of strong electronic 
correlations. These compounds include many cuprates, 
nickelates, vanadates and manganites, and are charac- 
terized by important and unique properties. In most of 
them the proximity of low-dimensional antiferromagnetic 
(AF) phases are the key to understand these properties. 

At the same time, the concept of spin laddeitl, orig- 
inally introduced to explain the presence of a spin 
gap in (VO)2P207 and later in layered cuprates like 
Sr„_iCu„-|-i02n (Refs. |2|j^) became an important the- 
oretical tool to understand the behavior of strongly cor- 
related systems. The physics of the two-leg spin ladder 
is characterized by the existence of a singlet-triplet spin 
gap and an exponential decay of correlation functions. 
The ground state which can be thought to a good ap- 
proximation as a product of singlets living on the rungs 
is now well understood. However, the above mentioned 
cuprates, the important case of Sri4Cu2404i, as well as 
many other compounds like CaV205, actually contain 
layers of coupled two-leg ladders. These ladders are cou- 
pled by frustrated interactions in a trellis lattice which 
make its study with analytical or numerical techniques 
quite difficult. In principle, in the absence of frustration, 
a reduction of the gap as the interladder coupling in- 
creases is expected. Eventually, the systetri becomes gap- 
less at a quantum critical point (QCP)Bu and for larger 
coupling it behaves essentially as a two-dimensional (2D) 
spin-1/2 square antiferromagnet. Much less is knowji for 
the trellis lattice, although a Schwinger boson studyQ sug- 
gests the transition from a spin liquid to a possible spiral 
order as the interladder coupling IILC) increases. Quan- 
tum Monte Carlo (QMC) studiesD of this frustrated sys- 
tem are hampered by the "minus sign problem" and the 



possibility of using this powerful technique is severely re- 
duced. In this sense, one of the objectives of the present 
work is to study a model in which the frustrated AF inter- 
ladder couplings of the trellis lattice are replaced by much 
simpler ferromagnetic (FM) couplings in a square lattice. 
We expect that some of the physics of the frustrated syfe 
tem can be captured by this effective simplified model.Ll 
Besides, there are compounds which consist of FM cou- 
pled ladders like SrCu203. The results of the present 
work could be relevant to other FM coupled gaped sys- 
tems like the dimerized chains in (VO)2P207.eI Previous 
studies have compared the behavior of AF and FM frus- 
trated and nonfrustrated coupled gapless spin systems 
(spin chains). Q 

Alternatively, a renewed interest in coupled ladders 
comes from the high-Tc cuprate superconductors them- 
selves. A numberj-Ojf recent experiments, mainly neutron 
scattering studiedly, indicate the presence of incommen- 
surate spin correlations which in turn have been inter- 
preted as coming from the segregation of charge carri- 
ers into ID domain walls or "stripes" leaving the regions 
between them as undoped antiferromagnets. There are 
several theoretical scenarios that hajie pcedicted or that 
attempt to explain this stripe ordeJiJii^ but the origin 
of this order and its relation to superconductivity is still 
controversial. In particular, the problem of stripe for- 
mation in a 2D microscopic model like the t-J or Hub- 
bard models is extremely difficult to study with analyt- 
ical or numerical techniques. In principle, the inclusion 
of charge and spin degrees of freedom is essential for the 
understanding of this problem. However, it has been sug- 
gested that assuming the presence of a stripe structure 
it is very instructive to study its magnetic propeddpa. by 
using a model with spin degrees of freedom only.ollj In 
this simplified model, the AF insulating regions between 
the charged stripes are considered as n-leg isotropic lad- 
ders caupled by an effective interaction. In one such 
studyll3, following the initial picture from Ref. |o|, the 
insulating regions were considered as 3-leg ladders cou- 
pled by AF interactions. However, a numerical study of 
the 2D t-J modelE£l, as well as early studies of charge in- 
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homogeneities in Hubbard and t-J modelgl^l, indicate the 
formation of "bond-centered" stripes, i.e. doped two-leg 
ladders alternating with undoped two-legiadders. Cou- 
pled spin two-leg ladders were also studiedlJ but its rele- 
vance to the physics of Cu-0 planes is relative since they 
miss the essential ingredient that the magnetic order of 
the AF slices is 7r-phase shifted as emphasized in Ref. 
Hence, the second motivation for the present work comes 
from the assumption that this 7r-phase-shift between lad- 
ders can be modeled by taking a ferromagnetic coupling 
between them. 

In summary, the purpose of the present work is 
to study ferromagnetically coupled two-leg ladders and 
compare their behavior with the case of AF coupling. 
If this model is considered as an approximation of AF 
systems in the trellis lattice, the FM coupling is an effec- 
tive interaction coming from the frustrated interactions 
between ladders. If this model is considered as an ap- 
proach to the stripe phase of the cuprates, the FM in- 
teraction comes from a collective effect determined from 
the competition of charge and spin degrees of freedom. 
In both cases, the results of the present study lead to 
predictions which can be tested experimentally. We use 
essentially numerical techniques like QMC (world-line al- 
gorithm) which allows us to reach low enough tempera- 
tures so as to capture ground state properties, and exact 
diagonalization with the Lanczos algorithm (LD), com- 
plemented by the continued fraction formalism to com- 
pute dynamical properties. 



II. QUASI-ONE DIMENSIONAL STUDY. 

To gain insight about the effects of FM interladder cou- 
plings we start from the case of FM coupled AF dimers 
which are the simplest systems with a spin gap. We are 
thus led to ID or quasi-lD systems which are much eas- 
ier to study from the numerical point of view. Besides, 
systems with a random distribution of AF and FM cou- 
plings have received some theoretical attention and their 
possible piwsical realization in SrCuPti_p 
discussed.EZi The Hamiltonian is given by 



lipOe has been 



H — Hdii 



where: 



J Sa;l • Sa;2, 



'Winter — ^ ^ Jinter^a^h^i^j^ayi ' ^b\j ^ (1) 



where a is a dimer index and i — 1,2 labels the sites in 
a dimer. J = 1 for simplicity. Periodic boundary con- 
ditions are imposed in the longitudinal direction. There 
are several ways of coupling dimers. We consider here 
the simplest case of dimers forming a FM-AF alternating 
chain (Fig. |l|(a)). Another possibility is that of dimers 



forming a two-leg ladder with FM leg and AF rung in- 
teractions. rThis second case has already been studied 
numericall}^ but it is not relevant for the problems we 
wish to address. Besides, we consider the case of FM- 
coupled AF plaquettes instead of dimers in which case 
we have a two-leg ladder with FM-AF alternating inter- 
actions along the legs (Fig. |l|(b)). 
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(a) (b) (c) 

FIG. 1. (a) Coupled dimers forming an alternating chain; 
(b) coupled plaquettes forming an alternating ladder; (c) cou- 
pled ladders in a square lattice. Full lines (dashed) correspond 
to AF (FM) interactions. 

Our first concern in this section is the behavior of the 
spin gap starting from the situation of isolated dimers 
or plaquettes. Using exact diagonalization we computed 
the spin gap A by substracting the energies in the S = 
and S = 1 subspaces on finite clusters with up to 24 sites. 
The extrapolation to the bulk limit was done using the 
law Aoo + bexp{—L/Lo)/L. The final result is shown in 
Fig. ||. We notice that in the limit Jinter —oo we 




FIG. 2. Singlet-triplet spin gap in the bulk limit of FM 
coupled dimers (circles) and plaquettes (triangles) as a func- 
tion of A = — Jinter- The cross indicates the spin gap for the 
spin-1 chain A = 0.41 divided by 4 (from Ref. (19)). 
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recover the cases of a spin-1 chain for the coupled dimer 
system and a spin-1 ladder for the coupled plaquette one. 
It is easy to realize (by solving a two-dimer system and a 
two-site spin-1 system) that the gap for the spin-1 chain 
is four times larger than the gap obtained by the coupled 
dimer system when Jinter ^ ~oo and the gap for the 
spin-1 ladder is twice larger than the coupled plaquette 
system in this limit. Thus, in the former case we obtain a 
gap Acd x4 — 0.410, coiiicident with the value already re- 
ported in the literature.Ej For the spin-1 ladder we would 
obtain a gap A^p x 2 — 0.290 ± 0.008, smaller fean the 
gap for the S=l chain as predicted theoretically.E3 Qual- 
itatively, the important feature here is that the gap de- 
creases monotonically from the isolated dimers (or pla- 
quettes) case as Jinter —oo. In the case of FM coupled 
dimers, a monotonic behavior could be guessed from the 
fact that this system continuously evolves towards the 
valence-bond-solid picture of a spin-1 chain in that limit. 
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FIG. 3. Dynamical structure factor as a function of the 
frequency for (a) AF and (b) FM coupled dimers for several 
momentum increasing from Qx = tt/IO to tt in steps of n/10 
from top to bottom. These results were obtained for the 20 
site chain with LD. The same for (c) AF and (d) FM coupled 
plaquettes. In this case, q = {qx,<ly), where x is the longitu- 
dinal direction, are shown in steps oi = 7r/4. These results 
were obtained for the 2x8 cluster with LD. 



On the other hand, the behavior of the coupled plaque- 
ttes system is not obviously predictable. 

The second point we want to examine is the behavior of 
the excitations of these systems, in particular the = 1 
excitations as can be measured by neutron scattering ex- 
periments. For this purpose, using conventional Lanc- 
zos techniques with the continued fractions formalism, 
we have computed the dynamical structure function (zz- 
component) ^(q, a;)Lil which is shown in Fig S. Already 
for the simplest case of coupled dimers (Fig ^(a),(b)) one 
can see an interesting feature which we will observe also 
for the coupled ladders case in the next section. The 
position of the gap which is at g = tt for AF dimerized. 
chains shifts to 7r/2 for the case of FM coupled dimers. E2l 
A similar behavior is observed for coupled plaquettes 
(Fig||(c),(d)), where the lowest energy peak changes from 
ilxTly) = ('"■,1') to (7r/2,7r) {x is the longitudinal direc- 
tion) by switching from AF to FM interplaquette cou- 
plings. In both cases, there is a transfer of spectral weight 
from the original AF peak to the FM one. This behavior 
is independent of the absolute value of Jinter, except for 
finite size effects. 



III. COUPLED LADDERS. 

We have now arrived at the central part of this 
work. The Hamiltonian for the system we consider now 
(Fig. ^c)) is essentially the same as (^ which we rewrite 
here for clarity: 



dde 



where: 



"Hladder — J ^a;Li ' Sa;/^i+i + J Sa-i^i ■ Sa-2,i, 

a.l.i 

"Winter — Jinter ^ ^ ^a\2^i ' ^a+l;l^ij (2) 



where a stands now for a ladder index and 1 — 1,2 for 
the two legs in a ladder. We have taken the case of an 
isotropic ladder by simplicity. J is again taken as the 
unit of energy. 

We start with the study of the ground state energy of 
the system of both FM and AF coupled ladders. To this 
purpose we have performed standard QMC simulations 
for L X L lattices with L = 4,6, 8, 12, and 16. Periodic 
boundary conditions in both directions are considered. 
For each lattice and set of coupling constants, we took 
T/J = 0.125, 0.100 and 0.07, and the Trotter number M 
at each temperature such that the error due to the time 
discretization is comparable or smaller than the statis- 
tical error. Typically, M = 140 for T/J = 0.07. We 
made runs up to 10® MC steps for both thermalization 
and measurement. We computed the energy in the total 
Sz — and Sz = 1 subspaces {Eq and Ei respectively). 
For each set of coupling constants, we extrapolated the 
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corresponding energies per site cq and ei to the bulk 
limit using the law Coo + bexp{—L/Lo)/L'^ in the gaped 
region and Coc + b/L^ in the gapless case. We obtained 
very close values for eoo.o and eoo,i which gives an in- 
dication of the good quality of the fits. The results for 
the ground state energy for FM and AF ladder couplings 
are shown in Fig. ^. An interesting feature can be no- 
ticed: the energies for both signs of Jinter are degenerate 
within numerical errors for \Jinter\ < 0.3. This situa- 
tion corresponds to a physics governed mainly by singlet 
dimers on the ladder rungs and in this case the sign of 
the coupling between these relatively isolated dimers is 
irrelevant. On the other hand, for \ Jinter\ > 0.3, for AF 
interladder coupling it is possible that the singlets delo- 
calize from a single ladder and finally form a "resonant 
valence bond" state or that the singlet-triplet excitations 
be replaced by j gapL ess magnon excitations. Previous nu- 
merical studiesrir'i precisely locate at Jinter ~ 0.3 the 
position of the QCP at which the ladder-like spin liquid 
is replaced by a long range 2D-like AF order thus choos- 
ing the second possibility. The important point we want 
to suggest is that in both cases the energy of the system 
would be lower than for the case of FM ILC where the 
singlets on the ladders still persist. To illustrate this sce- 
nario we have computed on the 16 x 16 cluster the spin- 
spin correlation functions S{r) = (SqS^) for r — (1,0) 
(leg direction), (0,1) (rung), and (1,1), inside a ladder, 
and r = (0, 1) between two ladders. These correlations, 
normalized in such a way that 5^0) — 1, are shown in ab- 
solute value in the inset of Fig. |4| as a function of | Jinter \ ■ 
The differences between FM and AF ILC appear in the 
(0,1) (rung) correlations, which remain stronger in the 
former case, and most importantly, in the interladder 
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FIG. 4. Ground state energy per site in the bulk limit of 
AF (open circles) and FM (full circles) coupled ladders as a 
function of the absolute value of Jinter- Error bars are smaller 
than the symbol size. In the inset, the intraladder (1,0) (cir- 
cles), (0,1) (squares), (1,1) (triangles), and interladder (0,1) 
(diamonds) spin-spin correlations are shown as a function of 
I Jinter I for the 16 X 16 cluster. Open (full) symbols correspond 
to AF (FM) interladder couplings. 



(0,1) correlation which increases faster in the latter. Of 
course, this correlation is negative (positive) in the AF 
(FM) case. 

This indication of a difference between the two ILC 
cases can be traced to a more intimate level which would 
also provide experimentally measurable features. To this 
end, let us examine now the static structure factor S'(q) 
obtained by Fourier transforming the spin-spin correla- 
tions obtained by QMC at the lowest temperature at- 
tained. In the case of AF ILC the peak is at (tt, tt) in all 
the range from the isolated ladder, which corresponds to 
the gaped "quantum disordered" region, to the isotropic 
square lattice, but the extrapolation of its intensity to the 
bulk limit becomes nonzero only for Ji^j~> Jinter, cr, i.e. 
in the "renormalized classical" region lBi3 For FM ILC, 
as shown in Fig. |5|(a) for Jinter = —0.2, the situation 
is qualitatively different. The peak in S'(q) is now at 
(tt, 7r/2), a feature which is similar to the one seen in the 
simpler cases of coupled chains and plaquettes. This be- 
havior has been found for all clusters considered, and for 
all Jinter < except for finite size effects: the smaller 
\ Jinter] the larger the size needed to reach the bulk be- 
havior. This is illustrated for the 4x4 cluster. 
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FIG. 5. Static structure factor obtained by QMC for 
Jinter = —0.2. (a) For various cluster sizes. The clusters 
considered are 4x4 (circles), 8x8 (squares), and 16 x 16 (tri- 
angles). The curve with full circles corresponds to the 4x4 
cluster and Jinter ~ —0.3. (b) For the 20 x 20 cluster at 
T = 0.2, 0.3, 0.4, 0.5, 0.6, 0.8, 1.0, 1.5 and 2.0. The curves have 
been multiplied by 2.5 for clarity. 

A highly nontrivial behavior is found if the tempera- 
ture dependence is analyzed. In Fig. ||(b)) the evolution 
of the structure factor for the 20 x 20 cluster and Jinter — 
-0.2 is shown at T = 0.2, 0.3, 0.4, 0.5, 0.6, 0.8, 1.0, 1.5 and 
2.0. In the high temperature region (T > 0.8) the peak 
is located at (tt, tt). As T is lowered the peak starts to 
shift towards the zero temperature peak (tt, 7r/2) which is 
reached at T = 0.3. We found almost no variation with 
cluster size of these two crossover temperatures at this 
value of Jinter- Thc fact that only at a finite tempera- 
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ture the peak of the magnetic structure factor starts to 
be incommensurate is reminiscent to the one first found 
in Lai.6-£!;Ndo.4Sra;Cu04 {x « 0.8) where a charge-stripe 
order is developed at Tc — 65K foUowecl-by a spin-stripe 
order at a lower temperature Tg — 50KJi3 

An important quantity to compute is the bulk limit of 
the peak of the structure factor. This quantity is related 
with the behavior of spin-spin correlations at the maxi- 
mum distance along the ladder direction (x-axis), along 
the direction transversal to the ladders (y-axis) and at 
the maximum distance of the 2D cluster. In the case of 
the AF square lattice, this latter quantity is proportional 
in the bulk limit to the squared staggered magnetization 
and it should be equal in that lip4t to the static struc- 
ture factor at momentum (7r,7r).E3 The finite size scal- 
ing of S'(7r,7r/2) is shown in Fig. ^(a) for Jinter- = —0.2 
and —0.6. We have attempted extrapolations to the bulk 
limit using both exponential and power laws. Due to the 
fact that clusters with L = 4, 8, 16 and L = 6,12 be- 
long to two different sets (which is more noticeable for 
large values of \Jinter\)^ the extrapolation procedure is 
not very reliable. However, as shown in Fig ^(a), one 
can conclude that ^(Tr, 7r/2) is zero for Jmter = —0.2 and 
nonzero for Jmter = —0.6. The finite size behavior of the 
spin-spin correlation at the maximum distance along the 
ladder direction, Smax,x, which is the one with smaller er- 
rors in our simulations, is similar to the one for 5(7r, tt/2) 
and the extrapolated values are also zero (nonzero) for 

Jinter — 0-2 {Jinter — 0.6). 
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FIG. 6. (a) S{TT,n/2) (open symbols) and Smax,x (full 
symbols) as a function of 1/L for Jinter ~ —0.2 (circles) and 
—0.6 (squares), (b) Singlet-triplet spin gap in the bulk limit 
of FM (open circles) and AF (full circles) coupled ladders as 
a function of the absolute value of the interladder coupling 
constant. In the AF case, the value of Ji„ter at which A = 
is taken from Ref. 23. The lines are guides to the eye. 

This crossover in the behavior of S{'k,tt/2) as a func- 
tion of Jinter poscs US with the question of the existence 



of a point analogous to the QCP in the AF ILC case. 
In the limit of Jinter —oo the coupled ladder sys- 
tem becomes equivalent to a system of AF coupled spin-1 
chains, where a fpite coupling is necessary to change to 
a gapless regimccJ To answer this question, let us now 
examine the behavior of the singlct-triplct spin gap as 
a function of Jinter- Although this is not a convenient 
quantity to compute with QMC since it involves a differ- 
ence between absolute values of the energies and then for 
large clusters the error becomes comparable to its value, 
we could get an indication of the presence or absence 
of a gapless region. The gap was computed for finite 
clusters and then extrapolated to the bulk using the Isfff 
Aoo + hcxp{—L/Lo)/L (or a/L^ for the gapless caseES). 
The results are depicted in Fig. |^(b). For the AF case, it 
can be seen a rapid decrease of A as Jinter is inCiypased, 
confirming earlier predictions and calculations. BSEj The 
gap vanishes at the QCP, Jinter,cr ~ 0.3. For FM ILC 
we also obtain a monotonically decreasing behavior, sim- 
ilar to the one found in the previous section for coupled 
dimers and plaquettes. The gap seems larger to that of 
AF ILC but it could vanish at Jinter ~ —0.4 within error 
bars. The calculation of other quantities like correlation 
lengths and/or using more powerful techniques should be 
necessary to obtain a reliable estimation for Jinter, cr- 





I 


(a) 


i 

ii i 






Jin.er=0.4 










i 1 1 ! J^^V^(0,;t) 


— 












111 L\ ^--^^^ 




i 




W \ 


/' \}Sj^ " 























0.0 1.0 2.0 3.0 0.0 1.0 2.0 3.0 4.0 

co/J co/J 



FIG. 7. Dynamical structure factor obtained by LD on 
the 4x4 cluster vs. frequency at several momenta q for (a) 
AF and (b) FM interladder couplings. From bottom to top 
the values of q are (7r/2, 7r/2), (tt, tt), (tt, 7r/2), (tt, 0), (■7r/2, 0), 
(7r/2,7i-), (0,7r) and (0,7r/2). 

The final part of our study which can eventually lead to 
a deeper understanding of the excitations involved in this 
system is the analysis of the dynamical structure factor 
S'(q, uj) which has been done with LD as in the previous 
section. In this case, we have to limit ourselves to some- 
what smaller clusters but we hope that the qualitative 
features we found will survive in the bulk limit. Results 
obtained for the 4x4 cluster are shown in Fig. 0. For AF 
ILC (Fig. |^(a)) the peak in ^(q, w) is located at (tt, tt), as 
expected in the bulk limit for an AF order. When a FM 
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ILC is involved (Fig. 0(b)) it can be seen that consider- 
able spectral weight is transferred to the peak at (tt, 7r/2), 
which becomes also the lowest energy excitation. Results 
for the 6x4 cluster are quite similar and it is quite re- 
assuring that these results are consistent with the ones 
obtained with QMC and shown in Fig. |[ 



IV. CONCLUSIONS 

In summary, we have numerically obtained some exact 
(except for extrapolation procedures) results for ferro- 
magnetically coupled systems, in particular two-leg lad- 
ders. Our main results are embodied in Fig. |^ and Fig. |^, 
i.e. at zero temperature the peak of the structure factor 
is located at (7r,7r/2) and it corresponds to the lowest 
energy excitation. Fig. |6|(a) also suggests finite values 
of this peak of the structure factor and the spin-spin 
correlation at the maximum distance along the ladder 
axis in the bulk limit for strong Jinter- There also two 
crossover temperatures, a higher one at which the peak 
starts to shift away from (tt, tt) and a lower one at which 
the peak reaches its zero temperature position. Besides 
the intrinsic interest for the theoretical understanding of 
spin- 1/2 ladder systems, we will try to emphasize in this 
section their possible relevance for realistic compounds 
containing ladders. As mentioned in the introduction, we 
should co^idcr in the first place ladder compounds like 
SrCu203EZI, Sri4Cu2404i (which upon Ga-doping and un- 
der pressure becomes superconducting^j) and CaV2 05. 
In these compounds, the ladders are coupled forming a 
trellis lattice. In the former case the interladder cou- 
plings are actually ferromagnetic. In the others, due to 
the frustrated nature of the AF ILC one could speculate 
that to some extent they could be modeled effectively by 
FM couplings. For all these compounds, then we predict 
that neutron scattering experiments would show peaks at 
{<lx, ly) = (tt, 7r/2), where the x (y) axis is in the direction 
parallel (perpendicular) to the ladders. 

As also suggested in the introduction, our results could 
be related to the striped structure which dynamically ap- 
pears in the Cu-0 planes of high-Tc cuprates. In this case 
we are able to trace the origin of the neutron scattering 
peaks observed away from (tt, tt) to an effective ferromag- 
netic interaction between 7r-shifted insulating spin lad- 
ders. In fact, this behavior can be observed already for 
the simplest case of ferromagnetically coupled AF dimers 
as shown in Section ||. In this case, it is easy to verify 
on small chains by LD or QMC that the peak moves 
continuously from g = 7r/2 to tt as some of the FM cou- 
plings are replaced by AF ones. If this picture could 
translate to coupled ladders, then one would be lead to 
the conclusion that some of the spin two-leg ladders are 
TT-shifted while others are in phase in order to reproduce 
the experimentally observed incommensurate peaks. An- 
other feature we want to emphasize is the temperature 
evolution of the structure factor (Fig. |5|(b)): the peak 



at (tt, 7r/2) is reached at a finite temperature and there 
is a range of temperatures in which an incommensurate 
peak is present. As mentioned in the previous section, 
this is reminiscent of the order in which charge and the 
spin stripes appear in the cuprates as the temperature is 
decreased. The fact that the spin gap possibly remains 
finite for somewhat strong values of | Jinter \ is also inter- 
esting for the stripe scenario of cuprates although in this 
case our results for the bulk limit are affected by large 
error bars. Of course, the question of to what extent this 
model of FM coupled ladders could apply to this scenario 
should come of detailed comparison with more realistic 
models like 2D t-J or Hubbard models. 
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